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Abstract 

Two NASA microwave radiometers, the satellite-borne GPM (Global Precipitation 
Measurement) Microwave Imager (GMI) and the aircraft-borne CoSMIR (the Conical 
Scanning Millimeter-wave Imaging Radiometer), measure vertically- and horizontally- 
polarized microwaves emitted by cloud particles and the Earth below, providing unique 
information on ice crystal properties in clouds. Their data reveal that non-spherical ice 
crystals are common and they fall in a preferred horizontally aligned orientation in 
convective and optically thick clouds especially near cloud top. 

A bin (particle-size-resolving) microphysical model with an ice crystal shape 
representation is developed to simulate the evolution of ice crystal properties (i.e., size, 
shape and orientation), where the radiation effect on microphysics (REM) is taken into 
account. Since REM represents the effect of all (e.g., both infrared and solar) radiation on 
ice crystal temperature, it relies upon the ice crystal proprieties that determine how an ice 
crystal receives radiation. Definitely, REM is different from the radiative effects that 
cause sensitivity at the microwave frequencies in the GPM and CoSMIR observations. 

Model results show that horizontally-oriented ice crystals grow faster than vertically- 
oriented (or spherical) ones due to REM. When both horizontally- and vertically-oriented 
ice crystals coexist in an air parcel, the model results show that the former grow by vapor 
deposition whereas the latter shrink by sublimation and disappear eventually. These 
modeling results are supported by the GMI data and the CoSMIR observations from 
MC3E (Midlatitude Continental Convective Clouds Experiment) and OLYMPEX 
(Olympic Mountains Experiment) on the prevalence of horizontally-oriented ice crystals. 


Moreover, the REM-induced precipitation explains the CloudSat observations of rare thin 


52 clouds in the tropical mid-troposphere as well as the common diamond dust in the high 
53 latitudes. 
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1. Introduction 

Ice clouds in the upper troposphere play a crucial role in the Earth’s radiation balance 
due to their widespread and long lifetime (e.g., Webster and Stephens 1980, Hartmann et 
al. 1984, Harrison et al. 1993, Zeng et al. 2009). However, they are not represented well 
in current weather and climate models (e.g., Zhang et al. 2005, Klein et al. 2009, Jiang et 
al. 2012). Field campaign and laboratory observations show that the ice crystals growing 
from the vapor phase are usually non-spherical (Pruppacher and Klett 1997). Since ice 
crystal shape is important in radiation computation (e.g., Liu et al. 1998, Yang et al. 
2005), how to better model and observe ice crystal properties (1.e., size, shape and 
orientation) is imperative to parameterize non-spherical ice crystals in climate models. 

Since aircraft penetration observations of ice crystals are expensive and limited over 
few geographic regions, the remote sensing of ice crystal shape is becoming important 
(e.g., Sassen 1974, Intrieri et al. 2002, Noel and Chepfer 2010, Neely et al. 2013, 
Skofronick-Jackson et al. 2015). In 2014, the GPM (Global Precipitation Measurement) 
satellite was launched (Hou et al. 2014, Skofronick-Jackson et al. 2015), where GMI (the 
GPM Microwave Imager) measures the intensity of radiation (via brightness temperature) 
at high frequencies (e.g., 166 and 183 GHz) that were added to GMI specifically for 
observing ice, snow and graupel properties. GMI’s microwave polarization data (e.g., 
166V, 166H) can be used to analyze the distribution of ice crystal shape (Skofronick- 
Jackson et al. 2008, Emory et al. 2013, Gong and Wu 2017). To support the analysis of 
satellite data, CoSMIR (the Conical Scanning Millimeter-wave Imaging Radiometer) that 
functions similarly as GMI was flown on the NASA ER-2 and DC-8 aircrafts during the 


Midlatitude Continental Convective Clouds Experiment (MC3E) and Olympic Mountains 
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Experiment (OLYMPEX) for clouds and precipitation observations (Jensen et al. 2016, 
Houze Jr. et al. 2017). In this study, the polarization data of GMI and CoSMIR are 
analyzed for the distribution of ice crystal properties. To better understand the GMI and 
CoSMIR data on ice crystal properties, a bin microphysical model that explicitly 
considers crystal shape and the radiation effect on microphysics (REM) (e.g., Zeng 2008) 
is developed to replicate the evolution of ice crystals in shape and size, where REM is 
different from the radiative effects that cause sensitivity at the microwave frequencies in 
the GPM and CoSMIR observations. 

The study was introduced in the 2016 PMM (Precipitation Measurement Mission) 
meeting briefly (Zeng et al. 2016) and completely herein by five sections. In Sec. 2, REM 
is incorporated into a numerical model for ice crystal shape simulations. In Sec. 3, 
numerical simulations are carried out to show how REM impacts ice crystal size and 
shape. In Sec. 4, model results are compared with GMI and CoSMIR data for model 
evaluation. Sec. 5 gives a summary. Except for special illustrations, all symbols are 


defined in the Appendix. 


2. Representing the radiation effect on microphysics (REM) 
a. The REM processes 

REM consists of a series of processes. That is, radiation impacts the surface 
temperature of cloud particles that in turn affects the saturation water vapor pressure 
around the particles and eventually alters the diffusional growth rate of the particles (e.g., 
Hall and Pruppacher 1976, Roach 1976, Stephens 1983, Wu et al. 2000, Zeng 2008). 
REM exists in two common weather phenomena: dew and frost in unsaturated air (e.g., 


Koeppe and De Long 1958). It also exists in ice clouds as a precipitation mechanism 
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(Zeng 2008). This section shows how to incorporate REM into a bin model for ice crystal 
shape simulations, beginning with REM for a single ice crystal. 

For a single ice crystal of mass m, its microscopic radiative ratio 77 is defined as the 
ratio of the infrared flux incident on ice crystal surface to the radiative flux emitted by the 
ice crystal as a blackbody at its environmental air temperature (Zeng 2008). Hence, 7 
approximately represents the ratio of the inward radiative flux to the outward one at 
crystal surface. When 7 = 1, the ice crystal grows just as predicted by the classic theory 
of cloud physics without REM (Fuchs 1959, Rogers and Yau 1989). 

If the atmosphere were a black body, 7~1. In fact, the atmosphere has quasi- 
transparent regions within the radiative spectrum (or atmospheric windows) and thus is 
not a blackbody (e.g., Stephens 1983). Hence, 77 usually deviates from one. As a result, 
radiation alters the ice crystal temperature and subsequently the saturation water vapor 
pressure around the crystal. An analytical expression of the ice crystal growth rate is 
derived as a function of 77 (Zeng 2008). 

Figure | illustrates the saturation water vapor pressure F; against ice crystal size when 
7 > 1 and 7 < 1, respectively. The dependence of E; (or ice crystal temperature) on 
crystal size originates in the different sensitivities of the radiative and conductive energy 
fluxes at crystal surface to ice crystal size. For a spherical ice crystal of radius r, for 
example, the radiative fluxes at crystal surface are directly proportional to 7° (or crystal 
surface area), whereas the conductive energy flux is proportional to r (or crystal size) 
only. 

When 77 < 1, E; decreases with increasing crystal size. Hence, an ice crystal can grow 


in an unsaturated environment, reaching a size of millimeters (see Sec. 3a for modeling). 
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This growth mechanism resembles the formation of a dew-drop (or frost) in an 
unsaturated environment if the dew-drop and its supporter (e.g., wide grass leaf) are 
treated together for a heat balance. 

Consider REM for two adjacent ice crystals shown in Fig. 1. When 77 > 1 (red lines), 
the larger crystal has a higher temperature and subsequently larger saturation water vapor 
pressure around it. As a result, the larger crystal sublimates with the resulting water vapor 
being deposited onto the smaller crystal. In contrast, when 7 < | (green lines), the smaller 
crystal has a higher temperature and subsequently larger saturation water vapor pressure 
around it. Therefore, the smaller crystal sublimates with the resulting water vapor being 
deposited onto the larger particle, which serves as a precipitation mechanism (Zeng 
2008). 

Consider multiple ice crystals that coexist in an air parcel. They compete for water 
vapor available and REM alters their competition for water vapor. As a result, REM 
brings about a nonlinear evolution of ice crystal spectrum. Zeng (2008) developed a bin 
(or particle-size-resolving) model to simulate how REM impacts ice crystal size for a 
given ice crystal shape. In this paper, the model is extended to simulate how REM 


impacts both size and shapes of ice crystals. 


b. Representing ice crystal shapes 
Consider an air parcel with ice crystals that moves with vertical velocity w. Its 
tendency of the relative humidity with respect to ice Hi, based on the conservation of 


energy and water, is derived as (Rogers and Yau 1989, Zeng 2008) 


dH, 
dt 


dm, 
= QT) 9,7) )—* (1) 
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where T is air temperature, myx G=1, 2, ...) denotes the mass of differently sized ice 


crystals for a given crystal shape k, and 


1| Leg g 
T)= : 2 
2.7) ee £ @) 
RT RL. 
T)= le 3 
2)=p [2 ee @) 
All ice crystals are classified into five shapes (or k =1, ... , 5): (1) a sphere of radius r 


and surface area S = 4nr’ for graupel, (2) a prolate spheroid of semi-major and minor axis 


lengths a and b for column-like crystals, whose surface area is 


2 a —b? 
S = 2b) b+ ————— arcsin ———— ], (4) 


a 
Va’ —-b° a 
(3) an oblate spheroid of semi-major and minor axis lengths a and b for plate-like 


crystals, for which 


Saf Z puieai se) (5) 
a’ —b b 

(4) a prolate spheroid for bullet-shaped crystals, and (5) a sphere for rosette crystals that 

consist of bullets. All of the crystal shapes have their own formulas of ice crystal mass, 

terminal velocity, dimensional relationship and ventilation factors against crystal size. 

The formulas are empirically obtained from field observations (e.g., Auer and Veal 1970; 

Heymsfield 1972; Heymsfield and Iaquinta 2000; see Zeng 2008 for computational 


details). 


c. Computing 7 
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The growth rate of an ice crystal, or dm/dt in Eq. (1), is derived as a function of 7 


(Zeng 2008), which is rewritten, after some algebraic operations, as 


dm 4nC(H,-1) Sé,oT” c. F. 
= + ‘_jlj-en ~- (6) 
dt H,A+B  KF,f,(H,A+B)\R8,T €,0T 


where o is the Stefan-Boltzmann constant, C the stationary diffusion shape factor that 
accounts for the influence of crystal shape on the water vapor field around a crystal (for a 


spherical crystal of radius r, C =r), 


A= ty A, 1 (7) 
FL foKT\R,T 
RT 
—— (8) 
and 
3 
Sie phe eee eee (9) 
RT é0l )4nCKF, f, 


(see the Appendix for the other symbol definitions). 

The ratio 7 in (6) is computed as follows. Let F* and F’ represent the upward and 
downward fluxes of infrared radiation in the atmosphere, respectively, and define the 
macroscopic radiative ratio of an air parcel as (Zeng 2008) 


Fr+F- 
ah ey 10 
1. 20T* ( ) 
Once atmospheric variables are given, 7; is obtained easily by using a two-stream 
radiation package to calculate F* and F’ (e.g., Fu and Liou 1992, Chou et al. 1995). 
The ratio 77, is used to determine 7, which is illustrated with Fig. 2. Consider an ice 


crystal that receives the upward radiation (F*), downward radiation (F°) and horizontal 


radiation (~oT*), where oT* represents the blackbody irradiance at air temperature T. 
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Since the atmosphere is not a blackbody, the upward infrared flux F* and the downward 
one F’ usually deviate from the horizontal one (or o7*). Their contribution to an ice 
crystal is approximated with (F*+F°) times the maximum horizontal cross section area of 
the crystal, which depends on the shape and orientation of the crystal. 

A horizontally-oriented plate-like ice crystal is discussed first. It is mimicked by an 
oblate spheroid of semi-major and minor axis lengths a and b. Since its maximum 
horizontal cross section area is 7a’, the spheroid receives the upward and downward 
infrared fluxes za’F* and a’F’, respectively, whereas the rest surface area of the crystal 
(or S-27a’) receives the horizontal infrared flux oT* from the crystal environment. Hence, 
the total radiative flux received equals oT (S-22a7) + ma?(F*+F). Substituting (10) into 
the total radiative flux to eliminate F*+F” and then dividing the resulting expression with 


oT'S yields the expression of 77. That is, 


2 2 
a (n.-)), (11a) 


ale 


where S is calculated with Eq. (5). When b << a, 7% 7. 

A vertically-oriented plate-like crystal is still mimicked with an oblate spheroid of 
semi-major and minor axis lengths a and b, but its maximum horizontal cross section area 
becomes zab. Following the same procedure to (11a), 77 for a vertically-oriented plate- 


like ice crystal is obtained, i.e. 


(n.-1). (11b) 


Consider a bullet-shaped crystal similarly. A horizontally (or vertically) oriented 
bullet-shaped crystal is mimicked by a prolate spheroid of semi-major and minor axis 


lengths a and b that has a maximum horizontal cross section area of zab (or mb). Thus, 
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n-1==2 G1) (12b) 


are obtained for horizontally- and vertically-oriented bullet-shape crystals, respectively, 
where S is from Eq. (4). 


Furthermore, the connection of 77, to 7 for a spherical ice crystal is obtained similarly 


as 
poly Zi) (13) 
1 = ) 1. 4 


Equations (11)-(13) can be summarized into a unifying expression or 
n-l=a(n,—-1), (14) 


where the coefficient @ varies with crystal size, shape and orientation and falls between 
zero and one. For a vertically oriented ice crystal with extremely long axis (b << a), a=0 
and thus 7 = 1. For a spherical crystal, a = 0.5. For a horizontally oriented crystal with 
extremely long axis (b << a), a= 1 and thus 7= 7. 

Although 7, and 77 are connected via Eq. (14), they have different applications. Since 
the macroscopic ratio 7; is independent of crystal properties, it can be treated as a 
property (i.e., the asymmetry between horizontal and vertical irradiances) of an air parcel 
and therefore can be used as a variable in cloud or climate models. In contrast, the 
microscopic ratio 7 varies with ice crystal size, shape and orientation, and thus is suitable 
to discuss the growth rate of individual ice crystals (e.g., Fig. 1). Since (77, - 1) has the 
same sign as (7 - 1), 7, usually takes the place of 77 in discussing REM except for 


specification. 
3. Numerical simulations of ice crystal properties 
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The parcel model in Sec. 2b is used to simulate how REM impacts the ice crystals in 
an air parcel of 7, = 0.5, where 7, = 0.5 corresponds to a cloud position with F ~ O (or 
near cloud top) approximately. In the present paper, the effect of solar radiation on 
crystals (e.g., Stephens 1983, Zeng 2008) is not discussed explicitly (or F;=0 is set). 
Instead, it can be discussed implicitly by tuning the magnitude of 77, to take account of Fs. 
Since there is a net radiative cooling (or net 7, <1) near cloud top even in daytime (e.g., 
Liou 1980) and the GMI observations at 166 GHz showed that the solar effect is about 
one third of the infrared one (Gong et al. 2017), the present simulations with 7, <1 are 
carried out to understand REM qualitatively. 

Six numerical experiments are carried out to test the sensitivity of ice crystal 
properties to REM (see Table | for the experiment summary). Except for specification, 
all of the experiments are set up as follows. All model variables are represented in double 
precision. Their prognostic equations are integrated with 2048 bins, using a time step of 
0.2 s with a second-order accuracy. The initial relative humidity with respect to ice is set 
to 100%. The air temperature T = -30°C, pressure p = 300 hPa, and 77, = 0.5. In addition, 


the vertical velocity w = 0. 


a. Contribution of REM to snow formation 

Diamond dust (or clear-sky precipitation) is common in the arctic regions especially in 
winter (e.g., Liu et al. 2012, Bennartz et al. 2013, Shupe et al. 2013). It involves neither 
riming nor aggregation. Since it is “simple” and 77, < | in the arctic regions (see Sec. 4a), 
it is simulated here using the bin model to show how REM impacts the shape and 


orientation of ice crystals. 
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The default experiment Ph addresses the evolution of horizontally-oriented plate 
crystals. Its initial ice crystal spectrum, shown by the thick line in Fig. 3, has an ice 
crystal concentration of 0.5 cm™ and an ice water content (IWC) of 3.6x10* g m®. If 7, = 
1, the ice crystals would maintain their initial spectrum without snow formation. 
However, if 7, = 0.5, the situation changes. Figure 3 displays the modeled broadening of 
ice crystal spectrum with 7; = 0.5, showing that precipitating particles form in a few 
hours. The spectral discontinuity near 200 um coincides with the emergence of dendritic 
extensions at that size. 

The average (half) crystal size for a given shape k is defined as 

a, = YiaM j, [2M 
J J 
where Mjx represents the mass of ice crystals in a bin of j and k. Figure 4 displays the 
average crystal size, IWC and H; against time in Ph, showing that the modeled process 
resembles the frost formation. When the average crystal size increases with time, REM 
leads to the decrease in the average crystal temperature and the saturation water vapor 
pressure around the crystals. As a result, water vapor deposits on ice crystals and the 
environmental relative humidity (or Hi) decreases with time. 

An interesting phenomenon in Fig. 4 is the dramatic increase (or decrease) in ice 
crystal size (or relative humidity) from hour 5.5 to 6.5. This dramatic change originates in 
the crystal dimensional relationship as well as the ventilation factors for mass transfer 
and thermal diffusion. That is, the ratio of crystal length to thickness and the ventilation 
factors increase with ice crystal size (e.g., Auer and Veal 1970; Pruppacher and Klett 


1997), which brings about the dramatic increase of ice crystal size after the ice crystal 
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size 1s longer than 1 mm. The dramatic change explains the formation of diamond dust in 


the dry arctic regions (e.g., Shupe et al. 2013). 


b. Selection of ice crystal orientations 

Two experiments PhPv and ChCv are carried out to test the selection of ice crystal 
orientations due to REM. The PhPv simulates plate crystals with two orientations. It takes 
the same setup as Ph except that half of the crystals are horizontally oriented and the 
other half are vertically oriented. Figure 5 displays the modeled spectra of vertically- and 
horizontally-oriented ice crystals. The horizontally-oriented ice crystals, just like those in 
Ph, grow to form precipitating particles in a few hours. In contrast, the vertically-oriented 
crystals shrink due to sublimation and disappear eventually (Fig. 5a), because the 
vertically-oriented ones have higher temperature than the horizontally-oriented ones 
when 77, <1 (see Fig. 2). 

The ChCv experiment simulates column crystals with two orientations. It takes the 
same setup as PhPv except for column crystals. Its modeled spectra are displayed in Fig. 
6, showing that horizontally-oriented crystals grow faster than vertically-oriented ones. 
The comparison of Figs. 5 and 6 indicates that plate crystals grow faster than column 
ones. In spite of their difference in growth rate, both experiments show that REM favors 
horizontally-oriented ice crystals, which explains the observed prevalence of 
horizontally-oriented ice crystals (e.g., Neely et al. 2013, Shupe et al. 2013; see Sec. 4c 


for a simulation evaluation). 


c. Selection of ice crystal shapes 
Ice crystal habit (or shape) varies with temperature and supersaturation (e.g., Magono 


and Lee 1966, Murray et al. 2015, see Pruppacher and Klett 1997 for review). Since 
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REM impacts crystal temperature and supersaturation via crystal size and shape, REM 
can affect ice crystal shape by selecting fast-growing crystals to survive. In this 
subsection, two experiments PhS and PhCh are carried out to test how REM selects 
crystal shapes. 

Experiment PhS simulates spherical solid ice crystals (or graupels if they fall into the 
parcel) mixed with horizontally-oriented plate ones. It takes the same setup as PhPv 
except that spherical crystals take the place of the vertically-oriented plate ones. Just like 
those in PhPv (or Fig. 5), the horizontally-oriented plate crystals grow to form 
precipitation while the spherical ones disappear. Figure 7 displays the average radius of 
spherical ice crystals against time. The average radius increases with time first and then 
decreases dramatically after the plate crystals reach a size of ~ 200 um with the 
emergence of dendritic extensions. These results indicate that horizontally-oriented plate 
crystals are selected to survive by REM when they coexist with spherical ones in an air 
parcel, which partly explains that non-spherical ice crystals are common in cirrus clouds. 

Experiment PhCh simulates plate crystals mixed with column ones by the same 
number concentration, where all of the ice crystals are horizontally oriented. The 
experiment is set up as PhPv except that the horizontally-oriented column crystals take 
the place of the vertically-oriented plate ones. Just like those in PhPv, the horizontally- 
oriented plate crystals instead of column ones are selected to survive by REM, because 


plate crystals grow faster than column ones (figure omitted). 


d. Effect of vertical velocity 
Vertical velocity impacts parcel relative humidity that in turn affects ice crystal 


growth. In this section, Exp. PhPvW is carried out to test whether vertical velocity alters 
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the selection of crystal orientations. The experiment is set up as PhPv except for a vertical 
velocity of w = 1 cm s‘!. Figure 8 displays the modeled evolution of the horizontally- 
oriented plate crystals, and Fig. 9 displays the average crystal sizes of horizontally- and 
vertically-oriented crystals. The horizontally-oriented ice crystals, just like those in PhPv, 
grow to form precipitation in a few hours whereas the vertically-oriented ones disappear 
in 1.5 hour, showing that vertical velocity does not alter the selection of crystal 
orientations. 

Figure 9 also displays the time series of IWC and crystal size in both PhPvW and 
PhPv, showing that vertical velocity increases IWC significantly. Interestingly, the 
average sizes of ice crystals in PhPvW and PhPv are close to each other, which exhibits 
the insensitivity of ice crystal spectrum broadening to vertical velocity. These different 
sensitivities of IWC and crystal size to vertical velocity suggest that REM, instead of 
vertical velocity, is the major contributor to the crystal spectrum broadening. The 
sensitivities also suggest that REM can be used to overcome the modeling problem of too 


narrow crystal spectrum that exists in the current cirrus models (Starr and Quante 2002). 


4. Comparison to GMI and CoSMIR data 

Since horizontally-oriented ice crystals are favored in an environment of 77, < | (see 
Table 1 for the experiment summary), their geographical distribution depends on the 
distribution of 7, < 1. In this section, the geographical distribution of 77; is analyzed first 
and then compared to that of GMI and CoSMIR 166V and 166H brightness temperature 


differences as an evaluation of ice crystal model simulations. 


a. Meridional variation of n- 
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The radiative ratio 7, varies greatly with atmospheric temperature, water vapor, 
sea/land surface temperature and clouds. However, its distribution for clear sky is 
relatively simple and thus provides a first-order approximation to understand thin cirrus 
clouds (including diamond dust)!. 

The ratio 7, in clear sky is computed, using Eq. (10) and the Goddard radiation 
package with the data of mean atmospheric states in the low, middle and high latitudes 
(Chou ef al. 1995). Figure 10 displays the vertical profiles of 7, in winter and summer in 
the high latitudes, showing that 77, is less than one in both seasons. Moreover, 77, in winter 
is much smaller than that in summer, because the underlying surface temperature in 
winter is much lower than that in summer. For comparison, Fig. 11 displays the vertical 
profiles of 7, in the Tropics and mid-latitudes, showing that 7, < 1 in the lower and 
middle troposphere but 7, > 1 in the upper troposphere in the Tropics and middle 
latitudes. Generally speaking, 77, decreases with increasing latitude, from 7, > 1 in the 
upper troposphere in the Tropics to 77, < 1 in the Arctic. This meridional variation of 77, 1s 
clear in physics, because the underlying surface temperature and subsequently the upward 
radiative flux F* decrease with increasing latitude. 

Furthermore, clouds change the vertical profile of 7, significantly. For thick clouds, 
their lower part blocks the upwelling infrared radiation from the underlying surface and 


consequently 7, < 1 near and above cloud top, no matter how high the underlying surface 


' Suppose that there is only one ice crystal in clear sky. Although the crystal contributes little to the 
atmospheric radiation transfer, the atmospheric radiation impacts the crystal temperature effectively. 
Hence, the growth of the crystal can be understood with 7, in clear sky. Similarly, if a thin cloud 
contributes little to the atmospheric radiation transfer, its crystal growth can be understood approximately 


with 7, in clear sky. 
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temperature is (Zeng 2008). Since the model results show that horizontally-oriented ice 
crystals are favored radiatively with 7, < 1, it is predicted that the frequency of 
occurrence for horizontally-oriented ice crystals is high near the top of thick clouds and 
the frequency in thin clouds increases with latitude. Such prediction on ice crystal 


orientation is evaluated with observational data next. 


b. CoSMIR data 

Passive microwave (PMW) radiometers such as CoSMIR and GMI measure the 
upwelling radiances emitted by the atmosphere, clouds, precipitation and the Earth below 
(e.g., Skofronick-Jackson and Johnson 2011). Their data can be used to retrieve 
hydrometeor particle information and rain rates (e.g., Kummerow 1993, Olson et al. 
2001a, Liu 2004). Since ice crystals in the atmosphere scatter Earth’s upwelling 
microwave radiation especially at high frequencies (> 85 GHz), the amount of 
microwaves scattered and thus received by PMW radiometers depend on ice crystals aloft 
(Skofronick-Jackson et al. 2015). If the ice crystals are non-spherical (e.g., plate-like) and 
oriented with respect to the PMW viewing angle, the microwaves received by PMW 
radiometers become polarized (e.g., Roberti and Kummerow 1999, Adams et al. 2008). 
Hence, the difference in radiance between horizontally- and _ vertically-polarized 
channels, especially at high frequencies (e.g., 166 GHz), should be able to measure the 
general orientation of ice crystals. Generally speaking, when the GMI brightness 
temperature (7) at 166 GHz for vertically polarized microwaves is higher than that for 
horizontally polarized microwaves, there are horizontally-oriented ice crystals aloft (e.g., 


Roberti and Kummerow 1999, Adams et al. 2008, Defer et al. 2014). 
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CoSMIR and GMI are both PMW radiometers with high frequencies at 166 GHz and 
operating with conical scans (thus at an approximate 53.6 degree viewing angle to the 
hydrometeors in the cloud) (Wang et al. 2007, 2008; Draper et al. 2015). The CoSMIR is 
loaded on the NASA ER-2 aircraft to observe clouds downwards with a footprint of 
2.5x4.3 km?. Its data from MC3E are analyzed next, providing an analysis procedure for 
GMI data. 

MC3E was a field campaign in south-central Oklahoma, conducted during April 
through June 2011 and jointly led by the NASA GPM Mission and the U.S. Department 
of Energy’s ARM (Atmospheric Radiation Measurement) program (Jensen et al. 2016). 
In the campaign, CoSMIR and several radars [including the HIWRAP (high-altitude 
imaging wind and rain airborne profile) radar] were in the high-altitude ER-2 aircraft to 
measure upwelling microwaves and hydrometeors (Heymsfield et al. 2013). 

CoSMIR data from MC3E are analyzed here to show the ice crystal orientation in a 
mesoscale convective system (MCS). Figure 12 displays the horizontal cross-section of 
radar reflectivities at 10 cm wavelength (or 3 GHz frequency) at an altitude of 3 km at 
2234 UTC May 23 from the ground-based NPOL (NASA S-Band Dual-Polarimetric 
Radar) located at 36.5N and 97.2W. The figure also displays the ER-2 flight track from 
2223 to 2241 UTC, 23 May 2011. This flight period is chosen for analysis because its 
observations cover both deep convection and anvil clouds. 

The HIWRAP radar and CoSMIR observations in the flight period are analyzed. 
Figure 13 displays the vertical cross section of Ku-band (2.8 cm wavelength) reflectivity 
(dBZ) from the HIWRAP radar that scans downwards, showing that there is deep 


convection, stratiform clouds and anvil clouds during the flight period. The figure also 
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displays the difference in the brightness temperature 7, at 165.5 GHz between vertical 
and horizontal polarization (referred to here as polarization difference). Since CoSMIR 
scans at an off-nadir-angle of 53.6 degrees, a distance shift of 5.5 km is applied to the 
CoSMIR data so that the CoSMIR data in the figure correspond to the radar data below. 
Figure 13 shows that the brightness temperature 7, at 165.5 GHz for vertical 
polarization is much higher than that for horizontal polarizations over stratiform clouds 
and anvil clouds, whereas 7, for vertical polarization is slightly higher than that for 
horizontal one over deep convection (Heymsfield and Fulton 1994, Gong and Wu 2017). 
Such difference in 7, suggests that a large fraction of ice crystals in non-convective 
regions are horizontally oriented but only a small fraction of ones in deep convection are 
horizontally oriented. Indeed, in deep convection particles are experiencing updrafts that 
(1) cause them to rotate almost randomly with no preferred orientation via turbulence and 
(2) detrain them near cloud top so that REM has insufficient time to select horizontally- 


oriented crystals. 


c. Comparison of COSMIR and GMI data 

The satellite-borne GMI functions similarly as the aircraft-borne CoSMIR. Its 
footprint at 166 GHz is approximately 4.1x6.3 km? that is larger than that of CoSMIR 
(i.e., 2.5x4.3 km? in MC3E and 1.3x1.9 km? in OLYMPEX)’. In this subsection, the GMI 
and CoSMIR data are compared over the OLYMPEX region, showing that GMI and 
CoSMIR measured the polarization difference consistently. 


2 The footprint of CoSMIR depends on not only the CoSMIR off-nadir angle but also the height of 


CoSMIR (or aircraft). To better match the GMI off-nadir angle, CoSMIR on the DC-8 aircraft took an off- 


nadir angle of 49 degrees in OLYMPEX, bringing about a footprint of 1.3x1.9 km? for the conical scan. 
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OLYMPEX is a field campaign to observe precipitation over the Olympic Mountains 
near Seattle from November 2015 through February 2016 (Houze Jr. et al. 2017). 
Different from MC3E, OLYMPEX addressed cloud systems that were mostly stratiform 
and also emphasized GPM core overpasses, providing unique data to compare the 
observations of GMI and CoSMIR. 

In the campaign, CoSMIR was loaded on the NASA DC-8 aircraft to measure 
upwelling microwaves with nine channels, two of which are at 165.5 GHz with both 
vertical and horizontal polarization. The two cloud systems observed on December 3 and 
5 are selected for analysis when the GPM core satellite (with GMI) and DC-8 (with 
CoSMIR) flew over the campaign region at the same time. Specifically, CoSMIR 
scanned from 14:33 to 16:59 UTC 3 December, between 46.49 and 49.59°N and between 
121.89 and 126.26°W while GMI flew over the region at 15:22 UTC 3 December; and 
from 14:18 to 15:38 UTC 5 December, between 46.45 and 48.61°N and between 122.4 
and 125.65°W while GMI flew over the region at 15:13 UTC 5 December 2015. 

The data from CoSMIR and GMI over the CoSMIR scan regions are displayed in Fig. 
14 for comparison, showing that the polarization difference over stratiform clouds is 
positive and becomes largest at horizontal polarization 7, of ~200K and smallest at 
~270K. Such positive polarization difference is consistent with that over the stratiform 
clouds during MC3E. 

Figure 14 also shows that the scattering of GMI polarization difference is narrower 
than that of CoSMIR. Such difference between GMI and CoSMIR data is attributed to 
their difference in footprint (or scan volume). In other words, the GMI data at 166 GHz 


represent the average of Tp over a larger area than those of CoSMIR. Although the GMI 
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data are different from CoSMIR in scan volume, they still catch the same main 
characteristics of polarization difference as the CoSMIR data and hence are used to 


statistically analyze the polarization difference from GMI. 


d. Statistics of GMI data 

The statistical analysis of polarization difference is carried out using three-year data of 
GMI and CloudSat data (Stephens et al. 2002), which resembles the MC3E analysis in 
Sec. 4.b except that GMI and CloudSat replace CoSMIR and the HIWRAP radar, 
respectively. A GPM data product of 2B.CSATGPM compiles the coincident data of 
GPM and CloudSat while CloudSat flew over GPM data pixels. Its 3687 coincident 
datasets, which spans from 8 March 2014 to 4 December 2016, are used for statistical 
analysis. 

One of the datasets is analyzed first as an example. Figure 15 displays the vertical 
cross-section of CloudSat radar reflectivity and the GMI polarization difference AT) 
(=166V-166H) at 166 GHz along the same CloudSat track over the Amazon on 13 
November 2014. The figure suggests that the polarization difference is positive over thick 
clouds but close to zero over thin cirrus clouds in the Tropics. Note that the polarization 
difference is close to 10K in stratiform cloud regions where it is expected that the ice 
particles are not subject to the turbulence in convective regions and thus can remain in a 
relatively fixed orientation. This nearly 10K difference was also true for the CoSMIR 
data. 

The statistical analysis of GMI data supports the AZ; results obtained from the 
preceding cases. To show the sensitivity of AZ; to latitude, clouds and the underlying 


surface, atmospheric columns (or pixels) are classified into four thicknesses: deep- 
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convective, thick, thin and extremely thin clouds (including diamond dust and clear-sky 
cases) when the maximum CloudSat radar reflectivity above 3 km is larger than 20 dBZ, 
between 5 and 20 dBZ, between -20 and 5 dBZ and less than -20 dBZ, respectively. The 
thicknesses are further distinguished by their underlying surface of sea and land. Finally, 
the thicknesses are divided into three geographical zones or latitudes lower than 25°, 
between 25 and 50° and higher than 50°. 

The data of GMI and CloudSat are analyzed for the probability distribution function 
(PDF) of AT; at 166 GHz in the Tropics first. Figure 16 displays the PDFs of tropical AT, 
over four kinds of clouds. The PDFs over the thin (both continental and oceanic) clouds 
are close to those over the extremely thin clouds, indicating that the 166 GHz channel 
cannot sense the contribution of the thin clouds to AJ». In other words, AJ; over the thin 
clouds comes mainly from the underlying surface (Adams et al. 2008) and thus is treated 
as the background here for brevity. 

Figure 16 also shows that the AT; over the thick clouds is significantly larger than that 
over the thin clouds (or background distribution), exhibiting a positive contribution of 
thick clouds to AT». Specifically, the thick clouds bring about an additional AT; increase 
of ~6K in the Tropics. Such positive contribution of thick clouds to AT» indicates that 
thick clouds prefer a horizontally aligned orientation of ice crystals near cloud top, which 
is consistent with the distribution of 7, < 1 near cloud top (Zeng 2008). 

Deep convective clouds, as shown in Fig. 16, contribute positively to AT} too, but less 
than the thick clouds, which supports the case analyses in Fig. 13. Interestingly, the PDF 


of AT» over oceanic convective clouds exhibits two peaks at 1.5 and 5.5K, which can be 
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attributed to two physical processes: the upward transport of randomly oriented (via 
turbulence) ice crystals to cloud top and the REM near cloud top, respectively. 

To show the meridional variation of AJ), Fig. 17 displays the PDFs of AT> over the 
thick and the thin clouds (or background) in the low, middle and high latitudes. The 
figure clearly shows that thick clouds contribute positively to the 166 AT; in all of the 
latitudes. Even though the difference in PDF looks smaller in the middle and high 
latitudes than that in the low ones, the horizontally oriented ice crystals in thick clouds 
may contribute the same to AT} in the middle and high latitudes as those in the low ones, 
because the background AT; increases with latitude (Fig. 17) and thus a large amount of 
spherical (or randomly oriented) hydrometeors in the lower part of thick clouds 
depolarize the upwelling background (polarized) microwaves in the middle and high 
latitudes (Adams et al. 2008). In short, the positive contribution of thick clouds to AT; 
exists in all latitudes and has no obvious meridional tendency, which indicates that thick 
clouds prefer a horizontally aligned orientation of ice crystals near cloud top and is 


consistent with the distribution of 7, < 1 near cloud top (Zeng 2008). 


e. CALIPSO data 

Field campaign and other satellite observations complement the preceding evaluation 
of ice crystal properties modeling. The observations of ground-based CAPABL (the 
Cloud Aerosol Polarization and Backscatter Lidar) reveal a large frequency of 
horizontally-oriented ice crystals in arctic thin clouds (Neely et al. 2013, Shupe et al. 
2013). Furthermore, the lidar and radar data of the CloudSat and Cloud-Aerosol Lidar 
and Infrared Pathfinder (CALIPSO) satellites show that diamond dust is common in the 


Arctic, especially in later winter and early spring (Liu et al. 2012). These observations are 
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consistent with the geographical and seasonal distributions of 7, < 1 in the high latitudes 
(or Fig. 10). 

The lidar polarization data of CALIPSO reveal the global distribution of horizontally 
oriented ice crystals in optically-thin clouds, showing that many horizontally oriented ice 
crystals exist in warm ice clouds (>-35°C) but few in cold ice clouds (<-35°C; e.g., Noel 
and Chepfer 2010, Zhou et al. 2012). Such sensitivity of oriented crystals to temperature 
can be explained by the vertical distribution of 77; in Fig. 11. In the Tropics and mid- 
latitudes, 7, <1 at altitude below 10 km (or temperature above ~ -35°C) and 77, >1 at 
altitude above 10 km (or temperature below ~ -35°C), indicating that REM can bring 
about many horizontally oriented ice crystals in warm ice clouds but few in cold ice 
clouds. In addition, the fraction of ice clouds containing oriented crystals increases with 
increasing latitude (Noel and Chepfer 2012), which can be attributed to the decrease in 7, 


with latitude with the aid of REM. 


f. CloudSat data 

The CloudSat observations of tropical clouds are consistent with the vertical 
distribution of 7;. Figure 18 displays the vertical profiles of thin and thick cloud fractions 
over land and sea in the Tropics (or latitudes lower than 25°). The fractions of thin and 
thick clouds are obtained by counting the pixels with radar reflectivity between -20 and 5 
dBZ and between 5 and 20 dBZ first and then normalized by their total pixel numbers 
above 2 km, respectively, where 2 km is chosen to avoid surface contaminated data (e.g., 
Marchand et al. 2008). 

In Fig. 18, the vertical distribution of thick clouds (with radar reflectivity between 5 


and 20 dBZ) exhibits the trimodal characteristics of tropical convection (Johnson et al. 
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1999). Since the fraction of thick clouds can be used to measure the detrainment of 
convection approximately, the strong detrainments at ~7.5 km and below 2 km 
correspond to the modes of deep and shallow convection, respectively (e.g., Riehl and 
Malkus 1958), whereas the strong detrainment at ~4.5 km corresponds to the third mode 
of cumulus congestus (Johnson et al. 1999). 

The vertical distribution of thick clouds cannot explain the distribution of thin clouds 
(with radar reflectivity between -20 and 5 dBZ) especially the maximum thin cloud 
fraction near 11.5 km, although thick clouds work as a water source of thin clouds. 
Instead, the REM-induced precipitation can explain the vertical distribution of thin 
clouds. Since 77, < 1 below ~ 10 km (see Fig. 11), the precipitation works as a water sink 
of thin clouds there, which explains the small fraction of thin clouds below ~ 10 km. In 
contrast, 7, > 1 above 10 km and thus there is no water sink like the REM-induced 
precipitation there, which explains the largest thin cloud fraction near 11.5 km even 
though the detrainment rate there is not large. 

The thin cloud fraction decreases with decreasing height from ~ 10 to 5 km, which is 
attributed to the variation in REM-induced precipitation rate. For a given 7, (<1), the 
REM-induced precipitation forms faster at lower height (or higher temperature) because 
the timescale of REM decreases greatly with increasing temperature or decreasing height 
(Zeng 2008). 

The thin cloud fraction reaches its minimum at ~ 5 km because the air temperature is 
O°C there. REM induces precipitation effectively only for non-spherical ice crystals and 
no ice crystals exist in thin clouds below 5 km. Besides, REM cannot bring about 


raindrops at the expense of small droplets (Roach 1976; also see Fig. 7 for a simulation of 
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spherical crystals). Hence, the lack of efficient precipitation formation below 5 km 
explains the minimum of thin cloud fraction at ~5 km. 

In addition to the tropical observations, the satellite observations of arctic clouds are 
consistent with the distribution of 77, too. The lidar and radar data of the CloudSat and 
CALIPSO satellites reveal a relatively high frequency of diamond dust in later winter and 
early spring (Liu et al. 2012), which corresponds to the smallest magnitude of 7, in the 
Arctic in winter (see Fig. 10). They also reveal a PDF peak of cloud top height between 7 
and 9 km, which corresponds well to the maximum of 7, near 8.5 km (or 300 hPa) in Fig. 
10. In summary, the distribution of 7, agrees well with the GMI observations of ice 
crystal orientation and the CloudSat/CALIPSO observations, indicating that REM plays 


an important role in the water cycle of the atmosphere. 


5. Summary 

Although REM plays an important role in the water cycle of the atmosphere, it has not 
been incorporated into the current cloud and climate models. In this paper, its effects are 
modeled against the data from satellites (GPM and CloudSat) and field campaigns 
(MC3E and OLYMPEX), which are summarized below. 

e A bin microphysical model with REM is developed to simulate the evolution of 
ice crystal properties (i.e., size, shape and orientation). Its results show that 
horizontally-oriented ice crystals grow faster than vertically-oriented (or 
spherical) ones when 77, < 1. If ice crystals with different shapes and orientations 
coexist in an air parcel with 7, < 1, horizontally-oriented ones grow while 
vertically-oriented (or spherical) ones shrink and disappear eventually, bringing 


about a selection of horizontally-oriented ice crystals. 
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The three-year data of GMI microwave polarization at 166 GHz reveal that 
horizontally oriented ice crystals are common in optically thick clouds in all 
latitudes especially near cloud top, supporting the model prediction of ice crystal 
shape because 77 < 1 near thick-cloud top. 

The CloudSat data show that the fraction of thin clouds in the middle troposphere 
is much smaller than that in the upper troposphere in the Tropics. Such vertical 
distribution of thin clouds is consistent with the model prediction, because the 
REM-induced precipitation (or 77, < 1 below 10 km) exists in tropical thin clouds 
in the middle troposphere but not in the upper troposphere (or 7; > 1 above 10 
km). In addition, the different signs of 77,-1 below and above 10 km (or -35°C) 
explain the CALIPSO observations of many and few horizontally oriented ice 
crystals in warmer and colder thin ice clouds, respectively (e.g., Noel and Chepfer 
2010, Zhou et al. 2012). 

In the arctic regions, 7, < 1 from the troposphere to stratosphere. The bin model 
simulations show that REM can bring about diamond dust (or clear-sky 
precipitation) even in a dry environment. The simulation results explain the high 
frequency of diamond dust and the prevalence of horizontally-oriented ice crystals 
in the arctic regions (e.g., Liu et al. 2012, Bennartz et al. 2013, Neely et al. 2013, 
Shupe et al. 2013). Such consistency in diamond dust between model and 
observations extends the GMI results into the arctic regions. 

The bin model simulations show that REM, instead of vertical velocity, broadens 


the ice crystal spectrum in cirrus clouds, which suggests that REM can be used to 
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overcome the modeling problem of too narrow crystal spectrum that exists in the 
current cirrus models (Starr and Quante 2002). 
REM will have two practical applications: (1) incorporated into the retrieval algorithm of 
GMI high-frequencies data for better precipitation products and (2) into multi- 
dimensional cloud/climate models to better replicate clouds and precipitation in the water 


cycle of the atmosphere. 
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Appendix List of Symbols 
a/b : semi-major/minor axis length 
C: stationary diffusion shape factor 
Cp : specific heat of dry air 


D: diffusivity of water vapor 
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é : partial pressure of water vapor 

£; : saturation vapor pressure over ice 

Jn/fo : ventilation factor for mass transfer/thermal diffusion 

Fs : solar flux absorbed by an ice crystal 

F*/F’ : upward/downward infrared flux 

FF: factor for the kinetic effect in heat/water vapor diffusion 
g : acceleration due to gravity 

H=e/E; : relative humidity with respect to ice 

Hic : critical relative humidity, see (9) 

K : coefficient of thermal conductivity of dry air 

Ls : latent heat of sublimation 

m : mass of an ice crystal 

MUna) : mass of ice crystals with semi-major axis length shorter than a 
N(a) : number of ice crystals with semi-major axis length shorter than a 
RJR, : gas constant for dry air/water vapor 

p : air pressure 

r : radius of a sphere 

S : surface area of an ice crystal 

t: time 

T : (air) temperature 

Ty : brightness temperature 

w : vertical velocity 


& : bulk absorption efficiency of an ice crystal for blackbody radiation 
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é : relative bulk absorption efficiency of an ice crystal for infrared radiation 

o: the Stefan-Boltzmann constant 

7 : microscope radiative ratio of individual ice crystals, see (11)-(13) 

nz : Macroscope radiative ratio of an air parcel, see (10) 

Pal Py : density of air/water vapor 

At : time step for numerical integration 

AT> : difference in brightness temperature between vertically- and horizontally-polarized 
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831 Table 1 List of the Numerical Experiments 
Experiment | Crystal Shape (orientation) | w Notes 
Ph + Plate (horizontal) 0 Default experiment 
PhPv + Plate (horizontal) 0 
+ Plate (vertical) Selection of crystal orientation 
ChCv + Column (horizontal) 0 
+ Column (vertical) 
PhS + Plate (horizontal) 0 
+ Spherical solid crystals Selection of crystal shape 
PhCh + Plate (horizontal) 0 
+ Column (horizontal) 
PhPvW + Plate (horizontal) Icms! Effect of vertical velocity 
+ Plate (vertical) 
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Caption 


Figure 1 Schematic diagram of REM for two adjacent ice crystals. The saturation water 
vapor pressure around an ice crystal E; varies with crystal size, where the red and green 
lines denote E; with 77 >1 and 77 < 1, respectively. When 77 <1, water vapor pressure e is 
lower (or higher) than £; of the smaller (or larger) crystal. As a result, the smaller crystal 
sublimates with the resulting water vapor depositing on and therefore increasing the 


larger crystal. 


Figure 2 Schematic diagrams of horizontally (left) and vertically (right) oriented plate 
crystals that receive the upward (F*), downward (F°) and horizontal infrared radiation 
(~oT’). Since F* and F usually deviate from oT‘, the horizontally- and vertically- 
oriented crystals have different energy budgets and therefore different surface 


temperature. 


Figure 3 Evolution of the mass density dM(na)/dina versus the half crystal size a in 
Exp. Ph for the horizontally-oriented plate crystals at 7, = 0.5, T=-30°C and p=300 hpa. 
The thick line denotes the initial spectrum and the time interval between lines is 30 


minutes. 


Figure 4 Time series of the average crystal size (top), the ice water content (middle) and 


the relative humidity with respect to ice (bottom) in Exp. Ph. 


Figure 5 Evolution of the mass density dM(na)/dlna versus the half crystal size a in 
Exp. PhPv for the vertically- (top) and horizontally-oriented (bottom) plate crystals. The 
thick lines denote the initial spectra and the time interval between lines is 30 minutes. 


Figure 6 Same as Figure 5 except for Exp. ChCv on column crystals. 
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Figure 7 Average crystal size versus time in Exp. PhS. The thick and thin lines 
represent the average size of horizontally-oriented plate crystals and the average radius of 


spherical ice crystals, respectively. 


Figure 8 Evolution of the horizontally-oriented plate crystals in Exp. PhPvW with a 


vertical velocity of 1 cms‘! 


Figure 9 Time series of the average crystal size (top) and ice water content (bottom) in 
Exps. PhPv (solid) and PhPvW (dashed line). The thick and thin lines in the top panel 
represent the average sizes of horizontally- and vertically-oriented plate crystals, 


respectively. 


Figure 10 Vertical profiles of the radiative ratio 7, in clear sky in the high latitudes 


during wintertime (thick) and summertime (thin line). 


Figure 11 The radiative ratio 7, versus height in clear sky in the Tropics (red), middle 


latitudes (green) and the arctic regions during wintertime (black). 


Figure 12 The horizontal cross-section of radar reflectivities at 10 cm wavelength from 
NPOL at 3 km altitude and 2234 UTC, 23 May 2011 during MC3E. The red line shows 
the ER-2 flight track from 2223 to 2241 UTC 23 May 2011, during which the CoSMIR 


data are analyzed. 


Figure 13 The polarization difference AT» at 165.5 GHz between vertical and horizontal 
polarization obtained from the CoSMIR conical scan at nadir (top) and the reflectivities 
from the HIWRAP Ku-band (2.8 cm wavelength) measurements (bottom panel) along the 


flight track shown in Fig. 12 during MC3E. 
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Figure 14 Microwave polarization observations of the CoSMIR conical scan at 165.5 
GHz (red) and GMI at 166 GHz (blue) over the OLYMPEX region on December 3 (left) 
and 5 (right), 2015, where a symbol corresponds to a scan pixel. The three black lines, 
representing the polarization difference (PD) of 0, 10 and 20K, respectively, are used for 


scaling. 


Figure 15 The vertical cross-section of the CloudSat radar reflectivity (bottom) and the 
GMI polarization difference AJ, (=166V-166H) along the same CloudSat track (top) 
when both CloudSat and GPM flew over the Amazon. This dataset starts over 7.38S and 


63.35W at 18:05:45 and end over 7.6N and 66.55W at 18:06:52 UTC 13 November 2014. 


Figure 16 PDF of the GMI polarization difference AT, at 166 GHz over land (left) and 
sea (right) in the Tropics. The blue, black, green and red lines represent atmospheric 
columns (or pixels) with the maximum radar reflectivity above 20 dBZ (deep convective 
clouds), between 5 and 20 dBZ (thick clouds), between -20 and 5 dBZ (thin clouds), and 


below -20 dBZ (clear sky and extremely thin clouds), respectively. 


Figure 17 PDF of the GMI polarization difference AT; at 166 GHz over land (top) and 
sea (bottom) versus latitude. The thick lines represent columns with the maximum radar 
reflectivity between 5 and 20 dBZ and the thin ones with the maximum radar reflectivity 
between -20 and 5 dBZ (referred to here as background) in the low (left), middle (middle) 


and high latitudes (right). 


Figure 18 The fractions of thick clouds (left, with radar reflectivity between 5 and 20 


dBZ) and thin clouds (right, between -20 and 5 dBZ) observed by CloudSat over the 


42 


899 Tropics. Their values are normalized by their total cloud pixel numbers above 2 km, 


900 respectively. 
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Figure 1 Schematic diagram of REM for two adjacent ice crystals. The saturation water 
vapor pressure around an ice crystal E; varies with crystal size, where the red and green 
lines denote E; with 7 >1 and 77 < 1, respectively. When 77 <1, water vapor pressure e is 
lower (or higher) than E; of the smaller (or larger) crystal. As a result, the smaller crystal 
sublimates with the resulting water vapor depositing on and therefore increasing the 


larger crystal. 
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912 Figure 2 Schematic diagrams of horizontally (left) and vertically (right) oriented plate 
913 crystals that receive the upward (F*), downward (F°) and horizontal infrared radiation 
914 (~oT"). Since F* and F’ usually deviate from oT, the horizontally- and vertically- 
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oriented crystals have different energy budgets and therefore different surface 
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919 ‘Figure 3. Evolution of the mass density dM(na)/dina versus the half crystal size a in 
920 Exp. Ph for the horizontally-oriented plate crystals at 7, = 0.5, T=-30°C and p=300 hpa. 
921 The thick line denotes the initial spectrum and the time interval between lines is 30 


922 minutes. 
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927 Figure 4 Time series of the average crystal size (top), the ice water content (middle) and 


928 _ the relative humidity with respect to ice (bottom) in Exp. Ph. 
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943 +Figure 5 Evolution of the mass density dM(Ina)/dina versus the half crystal size a in 


944 Exp. PhPv for the vertically- (top) and horizontally-oriented (bottom) plate crystals. The 


945 thick lines denote the initial spectra and the time interval between lines is 30 minutes. 
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967 Figure 6 Same as Figure 5 except for Exp. ChCv on column crystals. 
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Figure 7 Average crystal size versus time in Exp. PhS. The thick and thin lines represent 


the average size of horizontally-oriented plate crystals and the average radius of spherical 


ice crystals, respectively. 
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Figure 8 Evolution of the horizontally-oriented plate crystals in Exp. PhPvW with a 


vertical velocity of 1 cms". 
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Figure 9 Time series of the average crystal size (top) and ice water content (bottom) in 


Exps. PhPv (solid) and PhPvW (dashed line). The thick and thin lines in the top panel 


represent the average sizes of horizontally- and vertically-oriented plate crystals, 


respectively. 
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1003 +~=Figure 11 The radiative ratio 77, versus height in clear sky in the Tropics (red), middle 


1004 latitudes (green) and the arctic regions during wintertime (black). 
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1012 Figure 12 The horizontal cross-section of radar reflectivities at 10 cm wavelength from 
1013) NPOL at 3 km altitude and 2234 UTC, 23 May 2011 during MC3E. The red line shows 
1014 =‘ the ER-2 flight track from 2223 to 2241 UTC 23 May 2011, during which the CoSMIR 


1015 ~~ data are analyzed. 
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1021 Figure 13 The polarization difference AT) at 165.5 GHz between vertical and horizontal 


1022 polarization obtained from the CoSMIR conical scan at nadir (top) and the reflectivities 
1023 ~=from the HIWRAP Ku-band (2.8 cm wavelength) measurements (bottom panel) along the 


1024 flight track shown in Fig. 12 during MC3E. 
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1030 ~=Figure 14 Microwave polarization observations of the CoSMIR conical scan at 165.5 
1031 GHz (red) and GMI at 166 GHz (blue) over the OLYMPEX region on December 3 (left) 
1032 ~=and 5 (right), 2015, where a symbol corresponds to a scan pixel. The three black lines, 
1033 _ representing the polarization difference (PD) of 0, 10 and 20K, respectively, are used for 


1034 _ scaling. 
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Figure 15 The vertical cross-section of the CloudSat radar reflectivity (bottom) and the 
GMI polarization difference AT, (=166V-166H) along the same CloudSat track (top) 
when both CloudSat and GPM flew over the Amazon. This dataset starts over 7.38S and 


63.35W at 18:05:45 and end over 7.6N and 66.55W at 18:06:52 UTC 13 November 2014. 
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Figure 16 PDF of the GMI polarization difference AT; at 166 GHz over land (left) and 
sea (right) in the Tropics. The blue, black, green and red lines represent atmospheric 
columns (or pixels) with the maximum radar reflectivity above 20 dBZ (deep convective 
clouds), between 5 and 20 dBZ (thick clouds), between -20 and 5 dBZ (thin clouds), and 


below -20 dBZ (clear sky and extremely thin clouds), respectively. 
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Figure 17 PDF of the GMI polarization difference AT; at 166 GHz over land (top) and 
sea (bottom) versus latitude. The thick lines represent columns with the maximum radar 
reflectivity between 5 and 20 dBZ and the thin ones with the maximum radar reflectivity 
between -20 and 5 dBZ (referred to here as background) in the low (left), middle (middle) 


and high latitudes (right). 
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Figure 18 The fractions of thick clouds (left, with radar reflectivity between 5 and 20 
dBZ) and thin clouds (right, between -20 and 5 dBZ) observed by CloudSat over the 
Tropics. Their values are normalized by their total cloud pixel numbers above 2 km, 


respectively. 
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